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<^ i Abstract 
» . 

J3 ' The nonlinear Hartree equation describes the macroscopic dynamics of initiaUy factorized A^- 

boson states, in the Umit of large N. In this paper we provide estimates on the rate of convergence 
of the microscopic quantum mechanical evolution towards the limiting Hartree dynamics. More 
precisely, we prove bounds on the difference between the one-particle density associated with the 
^ ' solution of the A^-body Schrodinger equation and the orthogonal projection onto the solution of 

, the Hartree equation. 

00 
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^ '. 1 Introduction 

^ ■ We consider an boson system described on the Hilbert space L^(M^^) (the subspace of 

! consisting of all functions symmetric with respect to arbitrary permutations of the N particles) by 
a mean field Hamiltonian of the form 



N , TV 

j = l Kj 



HN = y -/^., + ^Y.^{xi-xj). (1.1) 



We will specify later assumptions on the interaction potential V . Note the coupling constant 1/A^ in 
front of the potential energy which characterizes mean-field models; it makes sure that in the limit of 
large the potential and the kinetic energy are typically of the same order, and thus can compete 
to generate nontrivial effective equation for the macroscopic dynamics of the system. 



We consider a factorized initial wave function 

N 



Ll{W^^) B Vjv(x) = JJv3(xj) for some <f G H\R^) (1.2) 



=1 
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with normalization ||(^||^2(]g3) = 1 (so that ||?/'Ar||L2(]fj3JV) = 1) and we study its time-evolution tpN,t, 
given by the solution of the N body Schrodinger equation 

idt'4'N,t = H^TpN^t with initial data ipN,o = V'Af- (1-3) 

In (jl.2p and in what follows we use the notation x = (xi, . . . ,xn) S M^^. 

Clearly, because of the interaction among the particles, the factorization of the wave function is 
not preserved by the time evolution. However, due to the presence of the small constant in front 
of the potential energy in (jl.ip . we may expect the total potential experienced by each particle to be 
approximated, for large N, by an effective mean field potential, and thus that, in the limit N — > oo, 
the solution ipj^^t of p.3|) is still approximately (and in an appropriate sense) factorized. We may 
expect, in other words, that in an appropriate sense 

N 

■0Ar,t(x) ~ JJ(^t(xj) for large iV. (1.4) 
i=i 

If (jl.4p is indeed correct, it is easy to derive a self-consistent equation for the evolution of the one- 
particle wave function ipt. In fact, it follows from (jl.4p that the total potential experienced by a 
particle at x can be approximated by the convolution {V * \ipt\'^){x), and thus that the evolution of 
the one-particle wave function ipt is described by the nonlinear Hartree equation 

idtipt = -Aift + {V *\ipt\'^)ipt . (1.5) 



To understand in which sense (jl.4p holds true, we need to introduce marginal densities. The 
density matrix 7Ar^i = \'4'N,t){'^N,t\ associated with ipp^^t is defined as the orthogonal projection onto 
TpN,t (we use here Dirac's bracket notation; for f,g,h e L^{W^), \f){g\ : L^{W^) ^^(M"') is the 
operator defined by \f){g\ih) = {g, h)i2 /). The kernel of ^N,t is thus given by 

7Ar,i(x;x') = V'Ar,t(x)V'Ar,t(x'). 

(k) 

For k = 1, . . . ,N, we define then the A:-particle marginal density 7]y ^ associated with 'ipN,t by taking 

(fc) 

the partial trace of 'j^^t over the last N — k particles. In other words, we define jj^ ^ as a positive 
trace class operator on L^(M^'^) with kernel 

7^,l(Xfc; x';.) = J dXAT.fe jN,ti^k,^N-k;^k^^N-k) ■ (1-6) 

Since HV'A^.t IIl2{r3JV) = 1, we immediately obtain Tr 7^^^ = 1 for all > 1, k = 1, . . . ,N, and t E R. 

By the choice of the initial wave function (jl.2p . at time t = we have 7^q = \ip){ip\^^. It turns 
out that p.4p should be understood in terms of convergence of marginal densities. For a large class 
of interaction potentials V, for every fixed A; > 1, and t G M, one can in fact show that 

l^Nl^\Vt){vr'' asiV^oo (1.7) 

where ^pt is a solution of the nonlinear Hartree equation (jl.Sp . The convergence (jl.7p holds in the 
trace norm topology. In particular, p.7p implies that for arbitrary k and for an arbitrary bounded 
operator J^^^ on L'^{R^^), 
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as N ^ oo. The approximate identity ()1.4p can thus be be interpreted as follows: as long as we 
are interested in the expectation of observables depending non-trivially only on a fixed number of 
particles, the iV-body wave function tp^^t can be approximated by the A^-fold tensor product of the 
solution (f>t to the nonlinear Hartree equation (jl.Sp . 

The first rigorous proof of (|1.7p was obtained by Spohn in [llj, under the assumption of a bounded 
interaction potential V. The problem of proving (jl.7p becomes substantially more involved for 
singular potentials. In [7], Erdos and Yau extended Spohn's approach to obtain a rigorous derivation 
of the Hartree equation (jl.Sp for a Coulomb interaction V{x) = const/|x| (partial results for the 
Coulomb interaction were also obtained by Bardos, Golse, and Mauser in [2j). In [4], the Hartree 
equation with Coulomb interaction was derived for semirelativistic bosons; in the semirelativistic 
setting, the dispersion of the bosons only grows linearly in the momentum (for large momenta), and 
thus the control of the Coulomb singularity is more delicate. In O [5l [6], models described by the 
Hamiltonian 

N N 

Hn = Y1 -A,^ + iV ^ A^3^y(A^^(xi - X,)) with (3 G (0, 1] 

j=l i<j 

with an A^-dependent potential were considered (in the one-dimensional case, A^-dependent potentials 
were considered by Adami, Golse and Teta in [Ij). These models are used to describe systems 
of physical interest, such as Bose-Einstein condensates. Assuming the interaction to be positive 
{V{x) > for all x G M^) and sufficiently small, the main result was again a proof of the convergence 
(jl.7|) : this time, however, ipt is a solution of the cubic nonlinear Schrodinger equation (with local 
nonlinearity) 

■ r, A I i2 -.1 ^ bo if < /3 < 1 

idtft = -A(ft + cr\(pt\ <ft with cj = I g^^^ if /5 = 1 • ^^-^^ 

Here bo = J dxV{x) and ao is the scattering length of V. The emergence of the scattering length ao 
for P = 1 (for all other choices of < /? < 1 the coupling constant is given by bo, which is the first 
Born approximation to Svrao) is a consequence of the short scale correlation structure developed in 
solutions of the Schrodinger equation, which, in the case /3 = 1, is characterized by the same length 
scale 0{\/N) as the scale of the interaction potential. 

The results described above have been obtained by extensions of the approach introduced by 
Spohn in which was based on the study of the BBGKY hierarchy 

'N-V 



+ 



N 



- j ^Trfe+i [V{xj - Xk+i),-iN^t 

/ .7 = 1 



for the evolution of the marginal densities 7]\rj, k = 1,. . . ,N (here Trfc_|_i denotes the partial trace 
over the {k + l)-th particle; this hierarchy is equivalent to the Schrodinger equation (|1.3p for ipN,t)- 
Because of the compactness of the sequence 7jvi) ^ ^ k, the proof of (jl.7p reduces to two main 

steps. The first step consists in proving that an arbitrary family of limit points {'ylx>\}k>i satisfies 
the infinite hierarchy 

k k 

idtlt]t = j;[-A.^,7S] + ^TVfc+i [V{x, - x,+i),72j'^] • (1-10) 
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The second step is a proof of the uniqueness of the solution of (jl.lOp . Since the factorized family 
7oo\ — Wt){'ft\®^, with Lpt determined by (jl.Sp . is a solution of the infinite hierarchy (jl.lOp . these 
two steps are sufficient to obtain (|1.7p . 

Despite its many successes, this method has some limitations. The main one, from our point of 
view, is that, because of the use of abstract arguments related to the compactness of the sequence 
7!^^, this technique does not provide any information on the rate of convergence of to \ipt){^t\®^ ■ 

In some cases, instead of comparing the solution of (|1.9p with the solution of the infinite hierarchy 
(jl.lOp . it is also possible to expand it in a Duhamel series and to compare it directly with the 
corresponding expansion for the factorized densities \ift){^t\®^- This approach (see leads to 
bounds of the form 

Tr |^(;)-|^,)(^,r|<^ (1.11) 

for all sufficiently small times |t| < Iq. The restriction to small times is needed to guarantee the 
convergence of the Duhamel expansion of the solution to (II. 9p . Iterating the arguments used to 
obtain (jl.lip . one can derive bounds of the form 



Tr 



(fc) 
TNt 



< 



Nl 



which hold for all t € M, but deteriorate very fast in time and are therefore not effective and not 
very useful. Next theorem, which is the main result of this paper, provides much stronger bounds 
on the difference between the true quantum mechanical evolution of the marginal densities and their 
Hartree evolution; in particular it shows that for every fixed time t G M, the error is at most of the 
order 0{N'^/'^). 



Theorem 1.1. Suppose that there exists D > such that the operator inequality 

V^(x) < D (1-A:,A 



(1.12) 



holds true. Let 



for some € 



N 



witl^ 



V'Af(x) = JJy5(a;j), 
1. Denote by ilJN,t - 



(1.13) 



ipN the solution to the Schrodinger 



equation hl.3\) with initial data tpN,o = tpN, o,nd let 7^^^ be the one-particle density associated with 
i/^N,t- Then there exist constants C,K, depending only on the norm of ip and on the constant D 
on the r.h.s. of such that 



Tr 



7S -Wt){vt\ 



< 



c 



Ari/2 



Here ipt is the solution to the nonlinear Hartree equation 

idtift = -Aift + (y * \ipt\'^)(pt 

with initial data ipt=Q = <~p- 



(1.14) 



(1.15) 



'^In what follows, for a function / we will always denote by its l? norm, while, for an operator A, ||A|| will 
mean its I?" operator norm. 
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Remark 1.2. The assumption on the potential V means that the most singular potential we can 
handle is the Coulomb potential V{x) = k/\x\. Note that our theorem applies both to the attractive 
(k < 0) and the repulsive case (k > 0). In particular Theorem \1.1\ implies the result obtained by 
Erdds and Yau in fTf. 

Remark 1.3. Note that under the assumption \1.12]) on the interaction potential V , the nonlinear 
equation ( lj.i5|] is known to be globally well-posed in H'^i^). This follows from the conservation of 
the mass \\(p\\ and of the energy 

8{^) = [dx |V(p(x)|2 + 11 dxdy V{x - y)\^{x)\^\v{y)\^ 



2 „ 

and from the observation that there exist constants ci,C2 such that 

Siip) <ci\\ip\\jjr{l + \\ipf) and y\\jj,<C2{£{ip) + \\ipf + \\ipf). (1.16) 
Both bounds can be proven using that, by 

dyV{x-y)\^{y)\^ <e\\V^f + e-^\\^f 



I' 



for all e > 0, uniformly in x G M^. 

Remark 1.4. Instead of (T7^ we will prove that 

\h'N\-\^t)M\HS<j^e^' (1.17) 

where \\-\\hs denotes the Hilbert- Schmidt norm. Although in general the trace norm is bigger than the 
Hilbert- Schmidt norm, in this case they differ at most by a factor of tw^. In fact, since \ipf ){ipt\ is a 
rank one projection, the operator A = — \ ipt){ipt\ can only have one negative eigenvalue Xmg < 0. 
Since moreover 

Tr {l^N]t-Wt){Vt\)=^ 

it follows that the negative eigenvalue of A is equal, in absolute value, to the sum of all positive 
eigenvalues. The trace norm of A is equal, therefore, to 2\\n(.g\ = 2||A||, where \\A\\ denotes the 
operator norm of A. Since \\A\\ < \\A\\hs, we immediately obtain that Tr\A\ < 2\\A\\hs- 

Remark 1.5. The bound (TTT^ is not optimal. As mentioned above, for short times and bounded 
potentials, the quantity on the l.h.s. of |j. j^[ ) is known to be of the order 1/N . Nevertheless Theo- 
rem \l.l\ is the first estimate on the rate of convergence towards the mean-field limit which holds for 
all times and remains of the same order N~^/'^ for all fixed times. 

Remark 1.6. Although, in order to simplify the analysis, we only consider the rate of convergence 
of the one-particle density 'y^^\ to \ipt){ipt\, our method can also be used to prove bounds of the form 



iVl/2 



Tr 

for all j,t,N and for j- dependent constants C{j),K{j). 

^We would like to thank Robert Seiringer for pointing out this argument to us. 
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In this paper we avoid the use of the BBGKY hierarchy and instead revive an approach, intro- 
duced by Hepp in [9] and extended by Ginibre and Velo in [8] , to the study of a semiclassical limit of 
quantum many-boson system^. This approach is based on embedding the A^-body Schrodinger sys- 
tem into the second quantized Fock-space representation and on the use of coherent states as initial 
data. The use of the Fock-space representation is in particular dictated by the fact that coherent 
states do not have a fixed number of particles. 

The Hartree dynamics emerges as the main component of the evolution of coherent states in 
the mean field limit (or, in the language of O [8], in the semiclassical limit). The problem then 
reduces to the study of quantum fluctuations, described by an A^-dependent two-parameter unitary 
group ly(N{t;s), around the Hartree dynamics. In O [8], Hepp (for smooth interaction potentials) 
and Ginibre and Velo (for singular potentials) proved that, in the limit N — > oo, the fluctuation 
dynamics U^it', s) approaches a limiting evolution U{t; s). This important result shows the relevance 
of the Hartree dynamics in the mean held limit (at least in the case of coherent initial states). It 
does not prove, however, the convergence ()1.7p of the one-particle marginal density to the orthogonal 
projection onto the solution of the Hartree equation, nor does it imply convergence results for the 
evolution of factorized initial sates. The problem of convergence of marginals requires additional 
control on the growth of the numbei0 of fluctuations generated by the evolution U]\[{t;s). This 
analysis, which, technically, is the most difficult part of the present paper (see Proposition 13. 3p . is 
ne'tto. Another novel part of our work is the derivation of convergence towards Hartree dynamics for 
factorized initial sates from the corresponding statements for the evolution of coherent states. 

Although we are mainly concerned with the dynamics of factorized initial data, the result we 
obtain for coherent states (see Theorem 13. ip is of independent interest, especially because, in this 
case, our bound is optimal in its A^-dependence (for coherent states, we show that the error is at 
most of the order 1/A^ for every fixed time). 

The paper is organized as follows. In Section O we define the Fock space representation of the 
mean field system, introduce coherent states and review their main properties. In Section [3l we 
consider the evolution of a coherent state and we prove that, in this case, the rate of convergence to 
the mean field solution remains of the order 1/A^ for all fixed times. Finally, in Section HI we show 
how to use coherent states to obtain information on the dynamics of factorized states, and we prove 
Theorem 11.11 

2 Fock space representation 

We define the bosonic Fock space over L^(IR^,dx) as the Hilbert space 

^ = L'^{M.\ dx)^=" = C e L2(m3", dxi . . . dXrr) , 
n>0 n>l 

with the convention L^(]R^)'^='^ = C. Vectors in T are sequences „>o of n-particle wave 

functions ip^"^^ G L^(IR3"). The scalar product on T is defined by 

(V'l,V'2) = ^('(/'i"'',^2"^>L2(R3n) = ^^^^2°^ / d^l---da;nV'i"^(xi,...,X„)^^"^(xi,...,Xn). 

n>0 n>l'' 

^Mathematically, the semiclassical limit considered in [9] [8] is equivalent to the mean field limit considered in the 
present manuscript. 

^Fluctuations around the Hartree dynamics will be considered as particle excitations and thus it will be possible to 
compute their number. 

more precise discussion of the results of [SJIH], and of their relation with our work can be found at the end of 
Section [S] 
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An A'' particle state with wave function '0jv is described on T by the sequence {'0^"^}n>o where 
^/;(") = for all n ^ N and ^/i^^^ = tpN- The vector {l,0,0,...}G^is called the vacuum, and will 
be denoted by Q. 

On !F, we define the number of particles operator J\f, by (AA^)(") = Eigenvectors of J\f 

are vectors of the form {0, . . . , 0, ■0^'"^ j 0, . . . } with a fixed number of particles. For / G L^(M^) we 
also define the creation operator a*(/) and the annihilation operator a{f) on JT by 



^"^=1 (2.1) 

(a(/)V^)(") {xi,...,Xn) = V^ J dxW)i^^'"^'\x,Xi,...,Xn). 

The operators a*{f) and a(/) are unbounded, densely defined, closed operators. The creation oper- 
ator a*(f) is the adjoint of the annihilation operator a{f) (note that by definition a(/) is anti-linear 
in /), and they satisfy the canonical commutation relations 

[a{f),a*{9)] = (/,5)l2(m3), [a{f),a{g)] = [a* (f) , a* (g)] =0. (2.2) 

For every / G L^(M^), we introduce the self adjoint operator 

</>(/) =a*(/)+a(/). 

We will also make use of operator valued distributions a* and ax {x G M^), defined so that 

a*{f)= [ dxf{x)al 

\ (2.3) 

«(/) = / dxf{x)ax 

for every / G L^(M^). The canonical commutation relations assume the form 

[ax, a*y\ = - y) [ax, ay] = [a* , a*] = . 
The number of particle operator, expressed through the distributions aj;,a*, is given by 

M = J a*ax ■ 



The following lemma provides some useful bounds to control creation and annihilation operators 
in terms of the number of particle operator J\f. 

Lemma 2.1. Let f e L'^{R^). Then 

Hm < 11/11 iiATV^v^ii 

||a*(/)^||<||/||||(Ar+l)^/^VII (2-4) 
||</,(/)V'||<2||/|| II (AT +1)^2^,11 

Proof. The last inequality clearly follows from the first two. To prove the first bound we note that 
||a(/)^|| < I dx|/(x)|||a,V'|| < dx\f{x)\^y' (^J dx||G,V||'^'^' 



(2.5) 



11/11 l|A^'/V|| . 
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The second estimate fohows by the canonical commutation relations (j2.2p because 



HfW + WffUf < Wff {\W"^\\ + Uf) = ll/f II (AA + if' ^Pf . ^^'^^ 



□ 

Given xj; £ we define the one-particle density 7^^^ associated with -0 as the positive trace class 
operator on L^(R^) with kernel given by 

7?^(^;y) = ^^;;^(^'«:«:^'V')- (2.7) 

By definition, 7^^^ is a positive trace class operator on L^(M^) with Tr7^"^^ = 1. For every A^-particle 
state with wave function V'Af £ Lg(R^^) (described on !F by the sequence {0, 0, ... , ipN,0, 0, . . . }) it 



is simple to see that this definition is equivalent to the definition (jl.6p . 
We define the Hamiltonian Hn on by (WtvV')^"^ = wS^V^"\ with 

n 1 " 

j=l i<j 

Using the distributions ax,a*, TC^ can be rewritten as 

Hn = j dxV c,alV r,ax + ^ j dxdyV{x - y)ala*yayax . (2.8) 

By definition the Hamiltonian Tljsf leaves sectors of J-' with a fixed number of particles invariant. 
Moreover, it is clear that on the A^-particle sector, Tl]\f agrees with the Hamiltonian Hf^f (the subscript 
N in Ti^ is a reminder of the scaling factor in front of the potential energy). We will study 
the dynamics generated by the operator TCn- In particular we will consider the time evolution of 
coherent states, which we introduce next. 

For / E L^(M^), we define the Weyl-operator 

W{f) = exp {a* if) - a{f)) = exp f / dx {f{x)al - l{x)a^)] . (2.9) 



Then the coherent state il^{f ) G J- with one-particle wave function / is defined by 
Notice that 



^ vn! 

n>0 n>0 

where /®"' indicates the Fock-vector {0, . . . , 0, /®", 0, . . . }. This follows from 

exp(a*(/) - a(/)) = e'll^ll'/^ exp(a*(/)) exp(-a(/)) 

which is a consequence of the fact that the commutator [a(/),a*(/)] = ||/|p commutes with a{f) 
and a*{f). From Eq. (|2.1U|) we see that coherent states are superpositions of states with different 
number of particles (the probability of having n particles in '(/'(/) is given by e^H-^H ||/|p"/n!). 

In the following lemma we collect some important and well known properties of Weyl operators 
and coherent states. 
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Lemma 2.2. Let f,g £ 

i) The Weyl operator satisfy the relations 

W{f)W{g) = W^((7)Ty(/)e-2*^™</'9) = W [f + g)e-' ^'^'^^^3) . 

a) W{f) is a unitary operator and 

wifY = wif)-' = w{-f). 

Hi) We have 

W*{f)a:,W{f) = a, + f{x), and VF*(/Kiy(/) = a* + 7(x) . 

iv) From Hi) we see that coherent states are eigenvectors of annihilation operators 

aM) = fi^Mf) a{g)^{f) = {9, f) L^^if) ■ 

v) The expectation of the number of particles in the coherent state ipif) is given by that is 

{ip{f)Mi'{f)) 



Also the variance of the number of particles in ipif) is given by (the distribution of Af is 
Poisson), that is 

- m),M^{f)? -- 



vi) Coherent states are normalized but not orthogonal to each other. In fact 

mf),i^{9)) = e-Hll/P+ll9iP~2(/,s)) ^ \{i;{f),^{g))\ = e--2y-ar . 

3 Time evolution of coherent states 

Next we study the dynamics of coherent states with expected number of particles N in the hmit 
— > oo. We choose the initial data 

tpi^Nip) = W{^ip)n iov £ H\R^) with \\ip\\ = 1 (3.1) 

and we study its time evolution 'il^{N, t) = e^^'^'^^ip{^/N (p) with the Hamiltonian TIn defined in ()2.8p . 

Theorem 3.1. Suppose that there exists D > such that the operator inequality 

V\x) < D{1 - A,) (3.2) 

holds true. Let T^''^ be the one-particle marginal associated with 'ip{N,t) = e~*^^*VF(\/iV¥')^^ (as 
defined in \2. 7\ )). Then there exist constants C,K > (only depending on the H^-norm of (p and on 
the constant D appearing in ^3. 2\) ) such that 



Tr 



rS^?,-|v'*)(^*||<^e^* (3.3) 



for all t G R. 

Remark 3.2. The use of coherent states as initial data allows us to obtain the optimal rate of 
convergence 1/N for all fixed times (while for the evolution of factorized N -particle states we only 
get the rate ; see ^.14\ ))- 

Proof. The proof of Theorem 13.11 will occupy the remaining subsections of section [3j 
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3.1 Dynamics IAn of quantum fluctuations 

By (|2.7p , the kernel of F j^^^ is given by 

\/ N \ 

It was observed by Hepp in [9] (see also Eqs. (1.17)-(1.28) in 0) that 



+ \j <ixdyV{x-y) {ipt{x)Lpt{y)ala*y + ipt{x)ipt{y)a..,ay) 
+ j dxdy y(x - y) al {Mvhl + 



It follows from (IHIil) that 



(3.4) 



(3.5) 

= UN{s;t) axUN{t;s) 
where the unitary evolution UN{t; s) is determined by the equatiorl^l 

idtUN{t]s) = C]y{t)l{N{t; s) and Z^Ar(s;s) = 1 (3.6) 

with the generator 

CN{t)= j dxVxal.Vxax + j dx {y*\ipt?') {x)alax + j dxdyV{x - y)Tp^{x)ipt{y)a*ax 



(3.7) 



+ ^{n,UN{t;OralUN{t;0)n) (3i 

+ ^^^ {n,UN{t;OyaxUN{t;0)Q) . 
V N 



In order to produce another decaying factor l/yN in the last two term on the r.h.s. of the last 
equation, we compare the evolution Z^jv(t;0) with another evolution Z//jv(t;0) defined through the 
equation 

idtUN{t;s) =ZN{t)UN{t;s) with Un{s]s) = 1 (3.9) 



^Note that, explicitly, Wjv(f,s) = W*(^/iV(^^)e-''^"(*-">W^(^/iV(^!>,). 
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with the time-dependent generator 

^N{t) = j dxV^a*V^aa. + j Ax [V *\Lpt?) {x)ala^ + j dxdyV{x - y)Tp^{x)ipt{y)ala^ 

+ \ I da^dy V{x - y) {ipt{x)Lpt{y)ala*y + '^t{x)'^{y)a^ay) (3.10) 



+ ^ y dxdy V{x - y) alayttyttx ■ 



From ()3.8p we find 

= j^{n,i(Nit;oya*yaMN{t;0)n) 

' ^n,l(N{t]0)*a*y [l(N{t;0) -Z^7v(i;0)) fl) + (Q, [UN{t]G)* -UN{t;G)*) a*yUN{t]Q)n 



+ ( ( Vt,UN{t; 0)*a, 0) - Z^jv(t; 0) ) 1^ ) + ( 1^, 0)* - UN{t; 0)* ) a,^/^(t; 0)0 )) . 



(3.11) 

Here we used the fact that 

(^n,UN{t;0)*ayUN{t;0)n'^ = <|l7,Z?Ar(t; 0)*o* Z?Ar(t; 0)0^ =0. 

This follows from the observation that, although the evolution Un^I) does not preserve the number 
of particles, it preserves the parity (it commutes with (—1)-'^). Multiplying (j3.1ip with the kernel 
J{x,y) of a Hilbert-Schmidt operator J over L^(M^) and taking the trace, we obtain 

Tr j(T^^\-\ipt){ipt\ 
= Jj^ j dxdy J{x,y){ayUNit;0)O,,axUNit;0)n) 

+ dxdyJ{x,y)^t{x){ayUN{t;0)n, (z^^(t; 0) - Z^7v(t; 0)) 0) 

+ ^ / dxdy J{x,y)Mx){(l^N{t;0) -UNit;0)) n,a*yUN{t;0)n) 



/N 

+ 



j dxdyJ{x,y)Tpt{y){alUN{t;0)n, (z^jv(t; 0) - Z^jv(t; 0)) Q) 
+ dxdyJix,y)^tiy){(uN{t;0) -l{Nit;0)) n,aMNit;0)n) . 
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Hence 



Tr J r 



^(1) 

N,t 



1 

< — 

- N 



+ 



dxdy \J{x, 



1/2 



dx\\a^UN{t;0)^f 
Z^Ar(t;0) -Z^Ar(t;0)j J] j dx\ipt{x)\\\a{J{x,.))UNit;0)n\\ 
Z^Ar(t;0) -Z?Ar(t;0)) J] j dx\iptix)\\\a*{Jix,.))UN{t;0)n\\ 
UN{t;0)-l(N{t;0)^n\\ J dy\ipt{y)\\\a*{J{.,y))UN{t;0)n\\ 
UN{t;0)-liNit;0))n [ dy\ipt{y)\\\a{J{.,y)pNit;OM 



and therefore 



Tr J r 



.(1) 
N,t 



N 
2|| J||hs 



N 
2|| J||hs 



\\iUNit;0) -l(N{t;0)M\ \\iM + l)^/^UN{t;0)n\\ 
\\{UN{t- 0) - UN{t; 0))n\\ + l)^/^UN{t; 0)O|| . 



The proof of Theorem 13.11 now follows from Proposition 13.31 Lemma 13. 8| Lemma 13.91 and from the 
remark that the trace norm can be controlled, in this case, by twice the Hilbert-Schmidt norm (see 
Remark 3 after Theorem II. ip . □ 

Proposition 3.3. LetUN{t] s) he the unitary evolution defined in kS.b]) . Then there exists a constant 
K, and, for every j £ N, constants C{j),K{j) (depending only on \\^\\h^ c-i^d on the constant D 
appearing in 13.^) ) such that 



for all Ip £ J- , and for all t, s G M. 



(3.12) 



Remark 3.4. Proposition 13.31 states that the number of particles produced by the dynamics Ujy of 
quantum fluctuations is independent of N and grows in time with at most exponential rate. This 
N -independence plays an important role in our analysis. Its proof requires the introduction of yet 
another dynamics ti^^^\ whose generator looks very similar to Ci^{t) but contains a cutoff, in the 
cubic term, guaranteeing that the number of particles is smaller than a given M. 

Proof. We start by introducing a new unitary dynamics with time-dependent generator C^^^\t) 
similar to £Ar(i) but with a cutoff in the number of particles in the cubic term. 
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3.2 Truncated dynamics 

For a fixed M > (at the end we will choose M = N), we consider the time-dependent generator 
^N^Ht) = j dx Va.a*Va.aa. + j dx {V*\<pt\'^) (x) a*aj; + j dxdyV{x - y)Tpl{x)iptiy)a*yax 

+ \ j dxdyV{x-y) {ipt{x)(ft{y)alal+lfl{x)Tp^{y)a^ay) 

( r (3.13) 

+ -= y dxdyV{x - y) al ((^t(y)ayx(AA < M) + ^t{y)x{M < M)al) a. 



+ ^ y dxdy V{x-y) a^alaya^ 



and the corresponding time-evolution U^\t;s), defined by 



idtUi^'^ {t;s) = 6P m]^'^ {t- s) with U^^'^ {s;s) = l 



Step 1. in the proof of Proposition \3. 

Lemma 3.5. There exists a constant K (only depending on and on the constant D in 13. ^) ). 

such that, for all N, M £ N, £ J^, and t,s £R 

{Ul^'\t; s)i^,M^U^^'\t; s)^l^) < (V, (M + 1)^'V) exp (4^' K\t - s\{l + . (3.14) 

Proof of Lemma \3. 51 To prove ()3.14p we compute the time-derivative of the expectation of (AA+ I)-'. 
It suffices to consider the case s = 0. We find 

= {U^^^'\t; 0)V, [iC^^'\t), (AA + iy]ui;'\t; 0)^) 

= lm j dxdyVix - y) ipt{x)^t{y){U^M"\f^ 0)^-, (AA + mU^^'\t; 0)V) 

Im / dxdyV{x - y)lpt{y){ui^'\t; 0)V^, [alayx{Af < M)a,, {AT + l)^>if ^(t; 0)V^) 



Using the pull-through formulae a^M = {M + l)ax, a%J\f = {J\f — l)a*, we find 

[al, {M + 1)^] = ^ {-l)\M + Ifal, [a., {M + 1)^] = {M + if a.. 

k=0 ^ ^ fc=0 ^ ^ 

As a consequence, 

i-i 

( I ) (-1)^ (<(AA + l)^a; + (AA + ifaluy 

E 

fe=0 

i-1 / -x i-i 



:a;, {M + 1)^] = (A {-if (<(AA + l)^a; + {M+ ifala;] 
k=o ^ ^ 

E (aA^<«;(A^ + 2)^ + (AA+ l)l<a;(AA + 3)t) , 

k=o ^ ^ 

[a., (M + ly] = J2 (i) + = E (i) (-^ + i)^«x.aaI. 

k=0 ^ ^ A:=0 ^ ^ 
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Therefore 

d. 
dt 

i-i 



Im / dxdyV{x-y)^t{x)My) 



i-i 



(3.15) 



To control contributions from the first term we use bounds of the form 

dxdyV{x - y) ^t{x)^t{y) {u'^'\t; 0)i;, {Af + l)ta>;(AA + 3)-2U^^'\t; 0)i;) 

< J dx\ipt{x)\\\a^{Af + l)-2U^^'\t;0)i;\\\\a^V{x - .)^^^^ 

< const sup (I Vix - yf\My)\'y ||(AA + 3)^Z^^^')(t; 0)V^f 

<K||(AA + 3)"*^4^)(t;0)V^f . 
Here we used that, by (13. 2p . 

sup / dyV'^ix - y)|v7t(y)P < DytWjji < constL>||v?||^i < K (3.16) 

X J 

is bounded uniformly in t (as follows from p.l6p ). Similar estimates are applied to the term con- 

k k 

taining J\f'2a%a*y{M + 2)2 . 

On the other hand, to control contributions arising from the second integral on the r.h.s. of 
(I3.15p . we use estimates of the form 

dx (t; 0)V^, a>(y(x - .)v.Ox(AA < M)(Ar + l)la,AAlz^if )(t; O)^') 

< j dx\\a^{M + l)-2U)^"\t-Q)n \\a{V{x- .)^t)x{M <M)\\ \\a:,M'^u\^"\t- 

< Ml/2 ||^(^ _ _)^^|| ||^Mi^CA'/)(i. 0)^11 ^ l)t<')(t; 

X 

<™i/2||(AA+l)^Z^if^)(t;0)V^f. 



This implies 
d 



- {U'^''^ {t; 0) V, (AT + lyU^^'''^ {t; 0)^^^) 



dt 

< K{1 + ^M/iV) ^ (<')(t; 0)^, (AA + 3)'=Z^^^')(t; 0)V^) 

< 4^- K(l + TM/iV) (Z^J^^) (t ; 0) V^, (AA + l)^' (t ; 0) V) . 

From Gronwall Lemma, we find (j3.14p . □ 
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Step 2. of the proof of Proposition \3.3\ 

3.3 Weak bounds on the dynamics 

To compare the evolution hlN{t]s) with the cutoff evolution l/(^^\t;s), we first need some (very 
weak) a-priori bound on the growth of the number of particle with respect to UN^t; s). 

Lemma 3.6. For arbitrary t, s G R and ij) ^ J-, we have 

{i^MNit] s)MUN{t; s)*V> < 6(V', iM + N + 1)11^) . (3.17) 

Moreover, for every £ £N, there exists a constant C (i) such that 

{^MN{t; s)M^^UN{t] s)» < C{t) (M + Nf^^) (3.18) 
{i^MN{t; s)N^'+^UN{t; sfij) < C{i) iij, (M + Nf'+\N + l)V'> (3.19) 

for all s G M and ip £ J-'. 

Proof of Lemma \3.(A Eq. (|3.19p follows from (|3.18p . In fact, assuming (|3.18p to hold true, we have 



{^MN{t-^s)M^'+^UN{t-sr^l:) 

< ^{%b,UN{t;s)M^'+^UN{t;sri^) + ^{^MN{t;s)M^'UN{t;sr^) 

< Die) {i^,i^^ + Nf'+H^^ + '^)i^) 

for an appropriate constant D{t). 

To prove (j3.17p and (j3.18p we observe that, by ()3.5p . 

UUt]s)NUN{t]s) 



(3.20) 



(3.21) 



(Recall that (j){(p) = a*{ip) + a{ip) = f dx{ip{x)a* + ip{x)ax))- From Lemma \2A\ and Lemma we 
get 



{ij,U*Nit; s)MUNit; s)V> < 2(V', W*{VNips)iM + N + l)WiVNips)'il^) 

= 2{ij,{M + Vn<I){^s) + 2N + (3.22) 
< 6{^p,{M + N +1)^P) 

which shows (I3.17|] . To complete the proof of (I3.18p . we define 

Xt^s = iM- \/iVe^^^(*-"V(¥'t)e~'^'^^*"'^ + N) . 
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Then, using the notation ad^(i?) = [-B,^], it is simple to prove that there exists a constant C > 
such that 

< C{N + Mf and ad^^^(AA) < C(N + M) for ah m G N. (3.23) 

By induction it follows that, for every £ G N, there exist constants D{i), C{i) with 

{Af + NfX^-'^ < D{1) {AT + Nf^ and X^j < C{£) {AT + Nf^ . (3.24) 

In fact, for £ = 1 (j3.24p reduces to (|3.23p . Assuming (j3.24p to hold for all i < k, we can prove it for 
i = k hy noticing that 



x^iH^ + Nfxl;' < 2{Ar + N)x^;;-\Ar + N) + 2\[xl;\Ar]\ 

k~2 2 



< 2{M + iV)Xg-2(AA + N) + 4^ ad^;^™(^) 



m=0 
k-2 



^t-s 



(3.25) 



m 
t-s 



< 2{M + N)Xf^^'\j\f + N) + A^C Xl'^,{Af + iV)2x, 

m=0 

< D{k) (M + Nf^ 

for an appropriate constant D{k), and that, by (j3.23p and (j3.25p . 

< CX^^;\Af + MfX^;^ < CD{k){M + Nf^ = C{k){M + Nf^ 
In (j3.25p . we used the commutator expansion 



m=0 ^ ^ 



in the second line, the bound (|3.23p in the third line, and the induction assumption in the last line. 
From (f3:2T]) and ([SlM]) . we obtain that 



< C{£){W{^ips)^, {M + Nf^Wi^ifs)^) (3.26) 

= c{e){ip, (M + ViV0(^,) + 2iv)2V) . 

Analogously to p.24p . it is possible to prove that, for every ^ € N, there exists a constant C{£) with 

{M + VNcPi^,) + 2Nf' < C{l){N + Nf . 
Eq. (IXTHI) follows therefore from (lOHl) . □ 



Step 3. of the proof of Proposition \K 
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3.4 Comparison of the Un and W^^^ dynamics 

Lemma 3.7. For every j € N there exist constants C{j),K{j) (depending only on j, on W'^Wh'^ o.nd 
on the constant D in such that 

{Unit; s)^lj,M' {uN{t; s)-U^i'\t- s)) ip) 



{1 + ,/M/N) 

and 

it; s)^, (uNit; s) it; s)) ^) 



" M3il + ^/MjN) V "^"^ V / ;i I 



for all ip & J- and for all t, s G M. 



Proof of Lemma \3. 1\ To simplify the notation we consider the case s = and t > (but the other 
cases can be treated identicahy). To prove (|3.27|) . we expand the difference of the two evolutions: 

{UNit;mM' (uNit;0)-U'^'\t;0)) ^) 



{UNit;0)tP,Af^UNit;0) ( 1 - Z^7v(t; 0)*Z^^*'^(t; 0) ) V) 



N 



i I ds {UNit;0)^,M^UNit;s) i^CNis) - C);'\s))U'^''''is;0)iP) 
ds J dxdyVix — y) 



'N 

X {UNit; Q)i^M'UNit; s)al (^t(y)aj,x(AA > M) + ^tiy)xiM > M)al) adA^^\s; 0)V) 
' ds I dx{a,UNit;s)*MmNit;0)ip,aiVix - .)^t)xi^ > M)aJA^^'^\s;Q)il:) 



N 
i 



ds / dx{adANit;sfMmNit;^)ij,xi^ > M)a*(y(x - .)^t)aJA^i'^\s;Q)iP) . 

(3.29) 
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Hence 



< 



N 



+ 



' ds / dx\\aMt; s)*NmM{t; \\a{V{x - .)^t)aM^ >M + l)u]^'\s- 0)^/;|| 
ds [ dx\\aMN{t]srU^UN{t;0)ma^V{x - .)^t)a,x{^f > M)ul^'\s-,0)i;\\ 



< 



sup 



\\Vix-.)^t\\ ds dx\\a,UN{t;syM^UN{t;0)i^\\ 



+ 



N 



:sup||y(x- .)(^t|| / ds / dx\\aMN{t;s)*Ar^UN{t;0)tp\\ 

X Jo J 

X \\a,M'/\iM > M)U^^'\s;0)i;\\ 
'ds\\M'/^UN{t;srMmN{t;Om \\Mx{M > M)ul^'\s;0)^\\ 

where we used (j3.16p once again. From Lemma 13.61 we obtain 

||AAi/2^jv(t; sYAf^UNit; 0)i;f = {Af^UN{t; 0)i^,li{t; s)NUN{t; s)*Af^UN{t; 0)4^) 

< 6{Af^UN{t; 0)4^, {Ar + N + l)Af^UN{t; 0)i^) 

<c{j){i;,{^^ + Nf^+H^^ + m 

< C{j)N^^+^{i;, {AT + . 
Therefore, using the inequahty x{-^ > ^) ^ {M /M)"^^ , we obtain 

{UN{t; Q)il^,N' (uNit; 0) - ui^'^Ht; 0)) i;) 

Jo 

< ||(Ar+ l)i+Vll l^ds {Uf\s;0)^, 0)^)^/^ • 

Finally, from (j3.14p . we conclude that 



X \\a,M'/\iM >M + l)U^^'\s;0)i;\\ 



< 



C 



(3.30) 



< C{j){N/My\\{Ar + ir+Vf / ds exp {K{j) s (1 + ./M/N)) 



< 



To prove p.28p . we proceed similarly; analogously to (I3.29P we find 



{U\,''\t; 0)^, (UN{t; 0) - U\^'\t- 0) ^) 



ds j dx{aMN[t-,syMm^^'\t;d)i>,a[V{x - .)ipt)x{M > M)aM'i''^ 

ds j dx{aJAN{t-sfMm^^'\t;^)^,x{M > M)a*{V{x - .)ipt)aMii'\s;Om 
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and thus 



r. ft ' (3.31) 

< -y= / ds\\Ar'/^UN{t;srAr^ui}'\t;Om ||AAx(AA > M)ui^'\s;0)i;\\ . 
Again, applying p.lSp and (13.14p we find 

{U^^^'\t; 0)i;,Af^ {u^{t; 0) - ui,^'\t; 0)) < C^^^^^^±^^ exp {K{j) t (1 + y^)) . 



Step 4- Conclusion of the proof of Proposition HOI 

From ([3221), and (I3l4l) we obtain, choosing M = N, 

{UN{t;s)i^,J\f'UN{t;s)^P) = {UN{t;s)^,M' {UN{t;s) -U^^'\t-s))^) 

+ {{UN{t; s) - U^^'\t; AA^'<'^(t; s)^;) 
+ {U^^\t-s)il,,M'u]^'\t-s)^) 
<C(j)||(AA+iy+Vfe^^^'^'*"''- 



□ 



□ 



3.5 Approximate dynamics UN{t; s) 

We now consider the dynamics UN{t] s), defined in (j3.9p by 

idtUN{t] s) = CN{t)UN{t] s) with Un{s]s) = 1 

with the time-dependent generator 

CN{t) = j dxV:ra*V^a^ + y dx (x) a*a^ + y dxdy ^(x - y)^(x)v9t(2/)a*aj: 

+ ^ y dxdy V{x - y) {ipt{x)ipt{y)a*^al + lfl{x)lpi{y)a^ay) (3.32) 



Lemma 3.8. There exists a constant K > 0, only depending on \\^p\\hi and on the constant D 
appearing in 13.^) . such that 

{UN{t; 0)n,Ar^UN{t; 0)n) < e-^* . (3.33) 
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Proof. We compute the derivative 



= 2Im / dxdyy(x -y)99t(x)(^t(y)(Z?jv(t; 0)17, [a>*(AA+l)3]Z?^(t; 0)0) 



4Im y dxdyy(x-y)(^t(j;)(^t(y)(Z^jv(t;0)rj, 

{ala*y{M + 1)2 + (AA + l)a*a;(AA +1) + {N + lfala*y) UN{t; 0)0) 
4Im j dxdyVix - y)iptix)^tiy){l(Nit;0)n, 

{{AT - l)a*.a;(AA + 1) + (AA + l)a*a;(AA + 1) + (AA + l)a*a;(AA + 3)) UN{t; 0)0) 
4Im / dxdyV{x - y)^t{x)^t{y){UN{t; 0)0, (3(AA + l)alal{M + 1) - 44a;) Z7iv(t; 0)0) . 



Therefore 



-(Z^7v(t;0)O,(AA + l)3z^;v(t;0)O) 



12Im j dxipt{x){a,^{M +l)UN[t]d)9.,a*{V[x - .)ipt){J\f +l)UN{t]Q)^) 
- 16Im / dxipt(x){a.jAN{t; 0)0, a*(y(x - .)ipt)UN{t] 0)0) . 



Taking the absolute value, we find 



^^{UN{t; 0)0, {N + l)3z^7v(t; 0)0) 



<12y dx|v9t(2:)|||a^(AA + l)^/iv(t;0)O|| ||a*(y(x- .)'^t)(AA + l)Z^iv(t;0)O|| 

+ 16y" dx|(^t(x)|||a^Z?7v(i;0)O||||a*(y(x-.)¥'0^Jv(i;0)O|| 

< 28sup||y(x- .)95t|lll(-AA+ l)^/2Z?iv(t;0)O||2 

X 

< C {UN{t; 0)0, (AA + l)^z7^(t; 0)0) . 
Applying Gronwall Lemma, we obtain (|3.33p . 



□ 



3.6 Comparison of the and Wa? dynamics 

The final step in the proof of Theorem 13. II is the comparison of evolutions generated by Uj^ and U]\f. 

Lemma 3.9. Let the evolutions U]\f{t; s) and h(j\f{t; s) be defined as in i3. 6]) and i3.9\) . respectively. 
Then there exist constants C,K > 0, only depending on lly^H/fi and on the constant D in 113. such 
that 



(llN{t;0)-UN{t;0))n 



< 



c 



-,Kt 



N 



(3.34) 
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Proof. We write 

= Unit- 0) (l - UNit; 0yUN{t; 0)) n 

= — J ds J dxUN{t;s)a^(f){V{x — .)(pt)axUNis;0)Q . 



Hence 



iV Jo 



ds 



Next, we observe that 

J dxal(j){V{x - .)(ft)axl^N{s;0)n 

= J dydx{ayUN{s; 0)U, - .)(pt)ayal(l){V{x - .)ipt)ajAN{s\ 0)J^) 

= j dydx{ayUNis;0)O,, (t){V (y - .)ipt)alay(j){V{x - .)(pt)ajAN{s;^)^) 

+ j dx{ajAN{s; 0)J^, 0(y(ar - .)ipt)^{V{x - .)<pt)aMN{s; 0)n) 

= J dydx{ayUN{s; 0)n, (a*(^(y(y - .)^t) + V{y - xp^ix)) 

X {^{V{x - .)ipt)ay + V{x - y)ipt{y)) ajAN{s; 
+ j dx{ajAN{s; 0)J^, 0(F(ar - .)^pt)(l){V {x - .)^t)ajAN{s; 0)0) . 

Therefore, we have 

j dx al4>{V{x - .)Lpt)axUN{s;{))9. 

= J dydx{axajAN{s; 0)0, ^(^(y - .)ipt)(l){V {x - .)ipt)ayajAN{s; 0)0) 

+ j dydxV{x - y)Tp^{x){ayUN{s\ 0)0, (j){V{x - .)ipt)ayaxUN{s; 0)0) 

+ J dydxV{x - y)ipt{y){axayUN{s; 0)0, (f){V{y - .)(pt)ajAN{s; 0)0) 



(3.35) 



+ j dydxV{x - y)^(pt{x)(pt{y){ayUN{s; 0)0, axUN{s; 0)0) 

+ j dx{ajAN{s; 0)0, 0(y(ar - .)^t)(l>{V{x - .)<pt)aMNis; 0)0) . 
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It follows that 

dxa*4>{V{x — .)(ft)(''xi^N{s;0)^l 
< snp\\V{x- .)iptf I dydx\\{M + 2)^/^a,ayUN{s]0)nf 

X J 

+ sup \\Vix - .)ipt\\ [ dydx\Vix - y)\\ipti^)\\\iM + lY'^ayUN{s- 0)O|| \\ayajAN{s; 0)J7|| 

X J 

+ suv\\V{y- .)^t\\ I dydx|y(x-y)||^t(y)|||a,aj,Z?^(s;0)O||||(AA + l)i/2a^.Z^^(s;0)O|| 
y J 

+ j dydxV{x - yf\^t{x)\\My)\\\ayUN{s;m\\\\axUN{s;m\\ 
+ sup||y(x-.)^if j dx\\{N +lf'^ajAN{s;mf ■ 
Using (|3.16p . we obtain 

j dx al(l){V{x - .)Lpt)axUN{s\ 0)17 



<C dydx\\a.^ayM^'^UN{s;G)VLf 



+ C { I dydx\V{x - y)\^\^t{x)\^\\ayM^/^UN{s;0)n\\^ 

1/2 



1/2 



dxdy\\ayaxUNis;0)O,\\'' 



+ Cl dydx\\axayUNis;0)Q\\' 



dydx\V{x - y)f\^t{y)maxAr^/^UN{s; 0)n\\' 



1/2 



1/2 



dydxV{x - yf\ipt{x)\'^\\ayUN{s;OM' 



+ C j dx\\axM^I'^UN{s;Q)nf . 



From 



dydx\V{x-y)\^\^t{y)\^\\axn^ < (^sup J dyVix - yf\^tiy)\'j WM'^'^pf < C\\M'/^i;f 
we thus find 

dx a*0(y (x - .)ipt)aMN{s; o)n ^ < c\\ (AA + if^^UNit; o)n\\^ . 



Inserting the last bound in (13.350 and using the result of Lemma 13.81 we obtain ()3.34|) . □ 
This concludes the proof of Theorem 13. li 

3.7 Discussion 

As mentioned in the introduction, our approach to the study of the mean field limit of the A^-body 
Schrodinger dynamics mirrors that used by Hepp and Ginibre-Velo in [9l [8] in the study of the semi- 
classical limit of quantum many-boson systems. In the language of the mean field limit, the main 
result obtained by Hepp (for smooth potentials) and by Ginibre and Velo (for singular potentials) was 
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the convergence of the fluctuation dynamics Z^Ar(t; s) (defined in (j3.6p ) to a limiting A^- independent 
dynamics U{t; s) in the sense that 



s- lim UNit;s) =U{t;s) (3.36) 

for all fixed t and s. Here the limiting dynamics U{t;s) is defined by 

idtU{t; s) = C{t)U{t; s) with U{s; s) = 1 

and with generator 

C{t) = j dx V^a*Va:ax + j dx [V *\ipt\'^) (x)a*.aa; + j dxdyV{x - y)lp^{x)(pt{y)ala^ 



+ 7; I 'ixdyV{x-y){ipt{x)(pt{y)ala*y + ipt{x)<ptiy)axay) 



(3.37) 



The convergence p.36p does not give any information about the convergence of the one-particle 
marginal associated with the evolution of the coherent initial state, to the orthogonal projection 

\ipt){ipt\- The definition of the marginal density T^''^ involves unbounded creation and annihilation 
operators. This also explains why the derivation of the bound (|3.3p in Theorem 13.11 is in general 
more complicated than the proof of the convergence (j3.36p . The proof of (j3.36p requires control of 
the growth of the expectation of powers of the number of particle operator M only with respect to 
the limiting dynamics. To prove (j3.3p . on the other hand, we need to control the growth of the 
expectation of M with respect to the A^-dependent fluctuation dynamics Z^Ar(t; s). 



4 Time evolution of factorized states 

This section is devoted to the proof of Theorem ll.il The main idea in the proof is that we can write 
the factorized A^-particle state ipN = (whose evolution is considered in Theorem 11. ip as a linear 
combination of coherent states, whose dynamics can be studied using the results of Section [3l 

Proof of Theorem \l.l[ We start by writing ijjN = ov, more precisely, the sequence 

{0, 0, . . . , 0, Vtv, 0, 0, . . . } = ^^^^^17 G 

as a linear combination of coherent states. While it is always possible in principle our goal is to 
represent ip^ with the least number of coherent states. 

Lemma 4.1. We have the following representation, 
with the constant 



/■Mr ^n 

n = dN — e'^^W{e^'^VNip)n (4.1) 
Jo 27r 
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Proof. To prove the representation (j4.ip observe that, from (j2.10p and since \\ip\\ = 1, 



d0_ 
2^ 



{a*{^)y 



n 



n. 



(4.3) 



□ 



The kernel of the one-particle density 7^^^ associated with the solution of the Schrodinger equation 



-UHn 



-n 



is given by (see ([27 



N 







27r 



(4.4) 







27r 



{Wie-'"^ VNip)n, a* (t)a^ (t)T^(e~^^2 ^Nip)n) 



where we introduced the notation ax{t) = e^^^^a^e ^'^'^ 



27r 



27r 



iWjvt^ ^-iHNt Next, we expand 



X (a^(t) -e-*^VAf93t(x)) W{e-''^WNip)Qj 



dN ftiy) 



2n 



d9^ 



27r 



27r 



d6i 



27r 



+ 



d%Vt{x) ['^dh rde2 ,e,NW,iN-i) 



N 



2tt 



27r 



X (VF(e^*''VA^99)Sl, (a*(t) -e^^VAf(?4(y)) Ty(e-*''VA^(/?)1^ ) 



+ (^t(x)(^i(y) 



27r 



d6'i 
2^ 



27r 



d^ 
27r 



L{Af-i) je2(Af-i) 



(4.5) 



We introduce the notation 
fN{x) =d% 



27r 



27r 



X (T^(e^*''VA^(^)J^, (a^(t) - e"'^''^ N^t{x)]W{e-'^''^ Nip)^l) . 



(4.6) 
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Since 



dN 



2lT 



— ( 

2-K 



,ie(N-i) 



°° / /-St Aft 



ie{N-i-j) j^j/ 



n 



e-^/2iV(^-i)/2 (a*(^))^-i 
dj\i , — — — i 1 



(4.7) 



VA^- 1! 



iV- 1! 



i>N-l 



we obtain, from (14. 5p . that 



N_ 

N 



2-K 



d6'i 
27r~ 



27r 



d6i 



2 -i9iN je2N 



2-K 



H \ 7^ \- Mxmiy) ■ 



Thus 



N 



< 



N 



2ix 



N 
2^ 



dS2 
2-K 



ay{t) - e-'^WN^tiy) ) W{e~''^WN^)n 



-102 



_^ \Mx)\\fN{y)\ ^ My)||/Af(x)| 



< 



r72 



27r 



d6'i 



27r 



27r 



aj,Z^^i(t; 0)1^11 ||a,Z^^2(t;0)Jl|| 



<^t(^)ll/^(y)l , \My)\\fN{x)\ 



N 



+ 



(4.9) 

where the unitary evolutions U^{t; s) are defined as in (13.61) . but with tpt replacecO by e~'^^ipt in the 
generator (j3.7p . Taking the square of (|4.9p and integrating over x,?/, we obtain 

dxdy |7j^J(x;y) - Lpt{x)i^t{y) 

<2^ /'^^^'^^ l|AAVX^(t;0)O|n|AAVX^(i;0)Of (4.10) 

+ dx|/^(x)|2 

Using Proposition 13.31 and the fact that c^at ~ A^^/^ to control the first term, and using Lemma 14.21 
to control the second term on the r.h.s. of the last equation, we find constants C, K, only depending 
on II^^IIj:^! and on the constant D in (I1.12|) such that 

C 



|HS 



< 



iVl/2 



exp(/Ct) . 



This proves (|1.17p and thus concludes the proof of Theorem II. 1[ 



(4.11) 



□ 



^We are making use here of the important fact that if (pt solves the nonlinear equation (|1.15p . then e'^ipt is also a 
solution of the same equation, for any fixed real 6. 
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Lemma 4.2. Let cpt be a solution to the Hartree equation lll.5\) with initial data G H^(M.^) with 
\\(p\\ = 1. Let 



X (W{e-''^WNip)n, (a^{t) - e-''^^VNipt{x)'j W{e-'^^VNip)n^ . 

Then there exist constants C,K (only depending on Hv^Hz/i CL^d on the constant D in iL12\) such 
that 

J dx\fN{x)f <Ce''K 

Proof. Using that 

^o^(t) - e-'^''VNipt{x)^ W{e-'^^^ip) = W{e-'^^^ip)U^^iO;t)a^U^^{t;0) 

where the unitary evolution U^{t;s) is defined as in (13. 6p . but with (ft replaced by e~'^^ipt in the 
generator (j3.7p . we can rewrite fN{x) as 



(V'(e2),<(0;t)aX'(^;0)^^) (4.12) 



with 

(■2-K 



^{92) = d% [ ^^e'''^''''^e-'''''W*{e-'''VNip)W{e~'''VN^)n. (4.13) 
Jo 27r 

Performing the integration over Oi, we immediately obtain 

^P{e2) = dN e-'^^^W*{e-'^^VNip)ip^^^-^^ . (4.14) 

It is also possible to expand ip{02) in a sum of factors living in the different sectors of the Fock space. 
From Eq. (|2.10p and Lemma 12.21 compute 



W*{e-''^^VNip)W{e-'^'VNip)n = W{-e-'^^y/Nip)W{e-''^^y/Nip)n 



-N Ne^''^-"!^ Sr^ \^ ^ ) ,„(S)m 

/ J I — r T 



(4.15) 



m>0 

which implies (using the periodicity in the variable ^i) 



m=0 

Switching to the complex variable z = e~*^^ we obtain 

m>0 



-\Tml2 r A-y 
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where the z integral is over the circle of radius one around the origin (in clock-wise sense). Changing 
variables z — > N z^ and using that djq = e^^'^^fW^./N^/'^ , we obtain 



M0^) = _ (iv - 1)! y il^e-ie2(m+i) / z-^e'(z - N^^p 



m=0 

oo 



m=0 

where we defined 



• 'mJ. 



(4.16) 



dz^ 

Comparing (j4.16p with (|4.14p . we obtain the identity 



'^m = ^^{e^{z-Nr)U=o- (4.17) 



y -Ml^ = 4 . (4.18) 



m=0 



It is also possible to obtain pointwise bounds on the coefficients TZm- From (|4.17p we deduce that 
for m<{N-l) 



m 



(' AT — 1 Vm W"~^ 777' 

The coefficients TZ^n turn out to be intimately connected with the classical system of orthogonal 
Laguerre polynomials. Recall that the associated Laguerre polynomial L^\x) admits the following 
representation 

" ^ ' ' k\{n - k)\{a + k)\ 

Therefore 

n^ = {-irm\L^^~^-^\N), 

which, for > m + 1, involves the value of the Laguerre polynomial L^\n) with a positive index 
a. Asymptotic expansions and estimates for the Laguerre polynomials is a classical subject, see |12j 
and references therein. However for the indices Q! = A^ — m — 1, n = m with N ^ m the value of 
X = N belongs to the oscillatory regime of the behavior of L^"^ (x) and the sharp estimates for those 
values of parameters have been only obtained recently in [10], where it is proven that, for a > — 1, 
n > 2 and the values of x E {q^, s^) the function lI^\x) obeys the bound 



.j-(a), M /(n + a)! x{s^ - q^) n 

where 

s = (n + a + 1) 2 + n2 , g = (n + a + 1) 2 — n2 , r{x) = {x — (f'){s'^ — x). 
As a consequence, we obtain that 
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Assuming that m < N and using the asymptotics {N — 1)! ^/^e ^ we obtain 

|L(^-™-i)(iV)|<m-i(m!)-^iVf 

and therefore 



< 



1 m 



_1 

m 4 . 



(m!)2 2 

bummarizmg, the coefficients Am = n„J{m^}/^N"'l^) appearing in the expansion (j4.16p of ^(^2) 
satisfy the bounds 

\Am\ < Cm~^/^ for all m < A and 

00 

Inserting (I4.16P into ()4.12p we obtain 



(4.20) 



fN{x) 



and therefore 



m=0 



(AA + 1)1/2 0)0 

d^ 



(4.22) 



27r 



(AA + 1)1/2 ^Y^(0;t)a^^/^(t; 0)17 



From ()4.20p . we obtain 



+ 1)3/2 E 1^™!'^ 



m=0 ' m=0 V ^ ^ - m>Ar 

On the other hand, from Proposition 13.31 we have 



const . 



(4.23) 



(AA + 1)1/2 Z^^(0;t)a^Z^^(t;0)J7 <^Ce^^ {M ^\faJU%{t;^)^ < Ce^* a^M'^U%[t;^)^ 



Thus, applying once more Proposition 13.31 we find 



2-K ^r\ 

dx |/,v(x)|2 < Ce^* / ^«(i; 0)17, AA5^/^(t; 0)0) < Ce^* . 



□ 
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